1 



Tail Asymptotics under Beta Random Scaling 



ENKELEJD HASHORVA0 and ANTHONY G. PAKEsi 

May 13, 2009 

Abstract: Let X,Y,B be three independent random variables such that X has the same distribution fimction 
' as YB. Assume that B is a beta random variable with positive parameters a, j3 and Y has distribution hmction 

■ H with H{0) = 0. Pakes and Navarro (2007) show under some mild conditions that the distribution function 
. ^a,/3 of X determines H . Based on that result we derive in this paper a recursive formula for calculation of 

H , if i?c(,/3 is known. Furthermore, we investigate the relation between the tail asymptotic behaviour of X and 
^ ' Y . We present three applications of our asymptotic results concerning the extremes of two random samples 

with underlying distribution functions H and Ha.p, respectively, and the conditional limiting distribution of 
CO . bivariate elliptical distributions. 

Key words and phrases: Beta random scaling; fractional integral; elliptical distribution; max-domain of at- 
' traction; asymptotics of sample maxima; conditional limiting results, estimation of conditional distribution; 

■ WeibuU-tail distribution; Gardes-Girard estimator. 

1 Introduction 

. Let A, y, B be three independent random variables such that 

> 

X = BY, (1.1) 

00 



(N 



where = stands for equality of the distribution functions. In our context the random variable B plays the 
role of a random scaling or multiplier. Clearly, if the distribution functions of Y and B are known, then the 
distribution function of A can be easily determined. In various theoretical and practical situations the question 
■ of interest is whether the distribution function of Y can be determined provided that those of A and B are 

known. Indeed, random scaling of y by i? is treated in several papers and different contexts, see for instance the 
recent contributions Tang and Tsitsiashvih (2003,2004), Jessen and Mikosch(2006), Tang (2006,2008), Pakes 
^ ' (2007), Pakes and Navarro (2007), Beutner and Kamps (2008a,b). 

b : 

^ Unless otherwise stated, in this article we fix B to be a beta random variable with positive parameters a, (3. 

If H denotes the distribution function of Y , then the distribution function of A (denoted by ffa,/?) is defined 
in terms of H and both parameters a, (3. If Y is another beta random variable, then A is the product of two 
such beta random variables, which have been studied extensively in the literature, see Galambos and Simonelli 
(2004), Nadarajah (2005), Nadarajah and Kotz (2005b, 2006), Dufresne (2007), Beutner and Kamps (2008a) 
and the references therein. 

Our main impetus for dealing with the beta random scaling comes from Pakes and Navarro (2007) which 
paves the way for the distributional and asymptotic considerations in this paper. Theorem 2.2 therein gives 
an explicit formula for the calculation of the distribution function H, provided that i?a,/3 satisfies some weak 
growth restrictions on its derivatives. Utilising the aforementioned theorem, we show in this paper that the 
distribution function H can be calculated iteratively without imposing any additional assumption on Ha^p. 
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This iterative inversion may lack the elegance of the explicit formula in Pakes and Navarro (2007), but it turns 
out to be quite useful in asymptotic contexts where we can define the tail behaviour of the survivor function of 
Y when that of the survivor function of X is known, and vice-versa. 

We present three applications of our asymptotic results: 

a) Determining which maximal domain of attraction contains H^jj when the membership of H is known; 

b) The derivation of conditional limiting results for bivariate elliptical random vectors; and 

c) New estimators for the conditional distribution function and the conditional quantile function of bivariate 
elliptical random vectors allowing one component of the random vector to grow to infinity. 

The paper is organized as follows. In the next section we give some preliminary results. The main result of 
Section 3 is the iterative inversion for -ffa,/3 ~ Theorem 13.31 below. In Section 4 we investigate the asymptotic 
relation of the survivor function of X and Y under conditions arising in extreme value theory, showing in 
particular that H is attracted to an extreme value distribution if and only if Ha.p is attracted to the same 
distribution. The direct implications are formulated (in Section 7) in a generality which subsumes the particular 
case of beta scaling. Conditional limiting results and estimation of conditional distribution function for bivariate 
elliptical random vectors is discussed in Sections 5 and 6. All proofs and some related results are relegated to 
Section 7. 



2 Preliminaries 



We introduce notation and then discuss some properties of the Wcyl fractional-order integral operator. A key 
result of Pakes and Navarro (2007) is recalled because it is crucial for our considerations. 

We use notation such as A" ~ to mean that X is a random variable with distribution function F, and 
F := 1 — F denotes the corresponding survivor function. The upper endpoint of the distribution function F is 
denoted by r^^ and its lower endpoint hy Ip- If a,/? > then beta(a,/3) and gamma(Q;,/3) denote respectively 
the beta and the gamma distributions with corresponding density functions 

(B(a,f3))-^x°'-\l-xf-\ xe(0,l), and -^a;"-iexp(-/3a;), x e {0,oo), 

r(a) 

where B{a,f3) is the beta function and T(a) is the gamma function. Since beta distributed random variables 
appear below in several instances, we use exclusively the notation i?a,/3 for a beta random variable with param- 
eters a, (3. On occasion it is convenient to extend the definition to understand P{i?o./3 = 0} = 1 if /3 > and 
P{Ba,i = l} = lifQ;>0. Unless otherwise stated, factors in products of random variables are assumed to be 
independent. 



Next, define the Wcyl fractional-order integral operator /3 > by 

1 f°° 

{lph){x) ■■= YX^J iy~xf-'hiy)dy, x > 0, (2.1) 
with h : [0,oo) — >iR a measurable function. The function /^ft, is well defined if (see Pakes and Navarro (2007)) 

x^^^\h{x)\ dx <oo 

is satisfied for all e > 0, in which case we write h €z with the understanding that (3 may assume negative 
values. Define further (consistently) loh := h. If /i is a density function of a positive random variable Y ^ H, 
then Iph is well-defined for every /3 > 0. Suppose 5 is a measurable function such that if Y ^ H, then 
E{Yf^-^\g{Y)\} < 00. Then we define 

{Jp^gH){x) = — - / iy~xf-'giy)dH{y), Va- G (Zj^, ), (2.2) 
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i-e., i7/3,g denotes the Weyl-Stieltjes fractional-order integral operator acting on the class of distribution functions 
oniR with weight function g. 

The Weyl fractional-order integral operator is closely related to beta random scaling. To see this, let a, /3 > 
and Y > and Ba,p be independent random variables such that 

X := YBaj:), where X ^ Ha,i3,Y ^ H, (2.3) 

and Ih > 0. In the light of equation (14) in Pakcs and Navarro (2007), for any x £ [In, th) we have 

H^A^) = r(^ + /?) x"(/;3P-o-^g)(x), (2.4) 
r(a) 

with ps the power function defined by 

Ps{x) ■= x'^ , s eM, X > 0. 

We mention in passing two important topics in probability theory and statistical applications where the Weyl 
fractional-order integral operator is encountered: a) the sized- or length biased law (see e.g., Pakes (2007), Pakes 
and Navarro (2007)); and b) the Wicksell problem (see e.g., Reiss and Thomas (2007)). For the essentials of 
fractional integrals and derivatives see Miller and Ross (1993). 

Now wc state three properties of/^/i. 

Lemma 2.1. Let /3,c be positive constants, and let h be a real measurable function. 

i) If h Cz 1f3+c, then 

Iplch = IJph = Ip+,h. (2.5) 

ii) Let 13" denote the n-fold derivative operator (n £lN). If the n-fold derivative ft.'-"-' := exists almost 
everywhere and h^"'^ G I/3, then 

D'^Iph = Iph^''^ (2.6) 

and 

D'^In = (-l)'7„_fc, k^l,...,n. (2.7) 
Hi) If X € (0,/?) and H is a distribution function onIR with H{0) = 0, then 

[I,3^xP~p{IxP~c.~xH)){x) = x-''{Ipp-o.-0H)[x), Vxe(0,oo). (2.8) 

The next theorem, which is an insignificant variation of Theorem 2.2 in Pakes and Navarro (2007) shows that 
the survivor function H can be retrieved by applying the differential and the Weyl fractional-order integral 
operator to Ha^p- 

Theorem 2.2. Let H jH^.p^OL^ (3 £ (0,oo) be as above, with iJa^^(G) = 0. If H^p'^^ is absolutely continuous 
and i?^"^ G Xg^a-ii i ~ 0, . . . ,n with 6 and n such that 

P + 6=:neN, (5 £ [0, 1), (2.9) 

then 

H{x) 



{-IT 



r(«) . 



r.a+0 



(IsDy^^Ha^Ai^) 



(2.10) 



holds for any x S (0, rn)- 
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3 Iterative Calculation of H 



Let X,Y, Ba^p, related by (|2.3p . be as above. In this section our main interest is the determination of H from 
the known form of Ha. p. As ah'eady mentioned, an cxphcit formula is presented as Theorem 2.2 in Pakes and 
Navarro (2007) (see ()2.10p above). If /3 G (0,1], then the only requirement for the validity of their theorem 
is that Ha^p{0) = 0, which obviously is fulfilled whenever H{Q) = 0. The following well-known multiplicative 
property of beta random variables is the key to our iterative version of Theorem 2.2 above. Specifically, if 
A e (0,/3), then 



Consequently, (|2.3p implies that 



Ba,P " Ba^\Ba+\,p~\- 



X = YBa^p = YBa^xBa+\,p~\- 



(3.1) 



Theorem 2.2 of Pakes and Navarro (2007) and (j3.ip implies the following result: 

Theorem 3.1. Let a, (3 he two positive constants, and let X, Y, Ba.p be independent random variables satisfying 
(HSl) with X - Ha,p, Y ^ H, and H{0) = 0. 
iJIfXe (0,/3), then 



HaA^) = r(« + /?) x"+^/^,AP-/3(/AP-a-Ag))(a;), Va;e(0,rH). 
r(a) 



(3.2) 



ii) // /3 — A e [0, 1), and 5 € [0, 1) is such that f3 — X + S = 1, then 



r(a + A) 



{a + X){Isp-a-x-iHa,p)i^) + iJ's,p_^_^Ha,p)i^) i Va;e(0,ri/). (3.3) 



We state next a simple corollary which is of some interest in the context of the Wcyl fractional-order integral 
operator. 

Corollary 3.2. Let H be a distribution function onlR such that H{0) ~ 0. Then for any x € (0,r//) we have 
x''-\jp,p_^_,^,H)[x) = {Dpal8P-a-i3H)ix) = -{DpaIpP-c.-pH){x). (3.4) 
Moreover, if H possesses the density function h, then 

{IpP-a-p+ih){x) = a{Lpp^a-pH){x) + x{LpD{p_a-pH)){x), xe{0,rH). (3.5) 



The main result of this section is the following iterative formula for computing H when Ha.p is known. 
Theorem 3.3. Let X ^ Ha,p, Y ^ H and Ba.p, a, (3 > be three independent random variables satisfying (|2.3p 
such that Ha,p{0) = 0. If (3q :~ (3 > f3i > ■ ■ ■ > f3k > Pk+i 0, with k S {0,1N} and Si,i < k + 1 are constants 
such that 



A, :=/3,_i-Ag (0,1], (5, :=1-A,, i = l,...,k + l, 
then we can construct distribution functions Hq := H, Hi, . . . , Hk+i = Ha,p such that 



(3.6) 



H^-l{x) 



r(a + A_i)' 



{a + l3,){Is^P~a~p^~iH,){x) + {Js^,p_^_,H,){x) , Vx e (0,rH).(3.7) 



Remark 3.4. (a) Let Bi^ Ba^^Pi^i >1 be independent beta random variables and independent of Y ^ H . If 
the random variable X with distribution function H„ ,n > 2 has the stochastic representation 



X = yY[b^\ c,e(0,oo), i-l,...,n. 



(3.8) 
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then Theorem \3.3\ implies that H can be retrieved recursively from i?„, provided that -ff„(0) ~ 0. 

(b) An interesting (open) question arises in connection with random products. Specifically, if M is a counting 
random variable taking positive integer values independent ofY,Bi,i > 1, such that 

X = where X^Hj^, (3.9) 

i=l 

then under what conditions on Af can we (recursively) compute the distribution function H if is known? 
Also arises a similar question if X, Y are related by 

X = Y[B3Bi + {1 - B3)B2]. (3.10) 



4 Tail Asymptotics 

The tail asymptotics of products have been studied in papers such as Berman (1983, 1992), Chne and Samorod- 
nitsky (1994), Tang and TsitsiashviU (2003, 2004), Jessen and Mikosch (2006), Tang (2006, 2008), and the 
references therein. Our asymptotic considerations below can be motivated by considering sample maxima. 

Specifically, let X;, K;, i = 1, . . . , n, be independent copies of X ~ YBa,p and Y , respectively, and 

Mx k ■= max X^, Myk max Fo, k>l 

l<j<k ' l<j<k 

be the corresponding sample maxima. From extreme value theory (see e.g., de Haan and Ferreira (2006), Falk 
et al. (2004, p. 23), Resnick (1987, p. 38)) if there are constants a„ > 0, 6„ such that 



lim sup 

ten 



H'^{ant + bn)-Q{t) = 0, (4.1) 



then we have the convergence in distribution 

{MY,n - bn)/an ^ My ^ Q , n ^ OO, (4.2) 

where Q is a univariate extreme value distribution (Gumbel, Frechet or WeibuU). If (14. ip holds (write H S 
MDA{Q)) it is of some interest to investigate the asymptotic behaviour of Mx,k,k > 1, where Xi,i < n are 
the results of a beta random scaling i.e., 

X, = Y,B,, B,^B^^p, i = 0,...,n, n > 1, (4.3) 

with Yi ^ H and the i?i'a and l^'s mutually independent. Thus Xi ~ H^^fs. A key question is whether Ha^p 
is in a maximal domain of attraction if H is, and conversely? We answer this below, as well as exposing the 
explicit tail asymptotic relations underlying (|4.3p . 

4.1 Gumbel Max-domain of Attraction 

If (|4.ip holds with Q = A the unit Gumbel distribution (A(a;) := exp(— exp(— x)), x GM), then there exists a 
positive measurable scaling function w (see e.g., de Haan and Ferreira (2006), Resnick (1987, p. 46)) such that 

H(x + t/w(x)) / N ^ /, ,N 

lim ^ = exp(-t), VteR (4.4) 

xlrn tl(x) 

is valid. We write H e AIDA{K, w) if (|4.4p holds. The scaling function w satisfies 

lim xw{x) = oo, and lim w{x)[rH — x) ^ oo, if rn < oo, (4.5) 
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and also the sclf-ncglccting property 



lim + = 1, (4.6) 

X^TH W{X) 



which holds locally uniformly for t G M; see e.g., Resnick (1987, p. 41). Note that most authors work with 
the so-called auxihary function l/w{x), but our convention follows Berman (1992) because results we prove are 
closely linked to some in his Chapter 12. 

Canonical examples of distribution functions in the Gumbel max-domain of attraction are the univariate Gaus- 
sian and the gamma distributions, which are special cases of distribution functions whose scaling functions have 
the form (for x large) 

1 + Li{x) 

where Li{x) is regularly varying at infinity with index 9fi,fi € (— oo,0) and r,6 arc positive constants. Note 
that 6* = 2 for the Gaussian case, and we have for the gamma(a, f3) case that 9—1, w(a;) = /3 and 

lim £i^±l) = exp(-/3i), yteM. (4.8) 

2;-^oo H{x) 



Distribution functions H that satisfy (|4.8p comprise what in other contexts is called the exponential tail class 
£(/3). See Pakes (2004) for references, and Pakes and Steutel (1997) where they are called medium-tailed. 

We state now the first result of this section, a close relative of Theorem 12.3.1 in Berman (1992); see Example 
1 below for the latter. In §7 we will state and prove the general proposition Theorem 1 7 . 41 which subsumes both 
direct assertions. 

Theorem 4.1. Let H^Ha^p he as in Theorem [XM Then H e MDA{K,w) iff (if and only if) Ha,p G 
MDA{A,w). If one of these holds, then 

HcAx) = {l+o{l))K{xw{x))-I^H{x), x]rH, (4.9) 

where K := V{a + (3)/T{a), and the density function ha.p of Ha.p satisfies 

lun^h-i(±_ ^ 1. (4.10) 

xUh w{x)H a,p{x) 

The asymptotic equivalence (|4.9|1 is the principal assertion here, as can be seen by noting that if one of the 
distribution functions F and H is in MDA{K, w) and they are related by 

F{x) = (1 + o{l))x'{w{x)YH{x), {x t rn), (4.11) 

where c, fi are real, then it follows from (|4.4p and (|4.6p that the other distribution function is in MDA{K, w). 

It is well-known that if iJ is a univariate distribution function with upper endpoint rn ^ oo and H G 
MDA{A,w), then H is rapidly varying (see Resnick (1987)) i.e., 

lim = 0, Vc> 1. (4.12) 

x^oo H{x) 

A necessary ingredient in the proof of Theorem 14.11 is the following rate of convergence refinement to (I4.12|) ; 
recall the first member of (|4.5p . 

Lemma 4.2. Let H he a univariate distribution function with rn = oo. If H G AIDA{K,w), then we have for 
any constant fi > 

lim {xw{x))'' = 0> Vol. (4.13) 
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Remark 4.3. (a) The self-neglecting property (|4.6p implies that the density function ha.p of Ha,p satisfies 

ha,fsix + t/w{x)) 

locally uniformly for t CzJR, provided that either H e il/_Dj4(A, w), or H^^js G MDA{K, w). 
(b) By Theorem \4-.l\ if Ha.p G MDA{A,w)^ then we can reverse (j4.9p obtaining 

H{x) = {l + o{l))-f^^{xw{x)fH^A^), xUh- (4.14) 

See Herman (1992) and Hashorva (2007d) for similar results. Further note that (j4.13p and (|4.14p imply for 
any c G (1, oo) that 

H{x) = o{H a^p{cx)), Ha^p{x) = o{H (cx)) , and {xw{x)Y H a^fj{x) = o(l), x t rn- 



We give next two illustrations of Theorem 17.41 

Example 1. (a) Theorem 12.3.1 in Berman (1992) follows from Theorem OTa) by taking (see (fTTTIl ) 

dp{u) = Pyi-Ba,,p>u} 

and checking that, since 1 — Ba^p = Bp^a, (|7.1ip holds with C = 2°' /oLB{a, (3) and the exponent (3 replaced with 
a. 

(b) Let 77, F be two distribution functions as in Theorem 14 . 1 1 and suppose that lu = Q and rn = oo. We assume 
that the random multiplier B has the stochastic representation 

B ^ XUi + (1 - X)U2, Ae(0,l), 

where Ui , U2 are two independent positive random variables such that for i = 1 , 2 

P{f/, > 1-s} = (l + o(l))c,s''% c„d^e(0,oo), siO. 

It follows that as s j 

PIS >i- si = (i + o(i)) r(di + i)r(d2 + 1) 

/-ii^>i ,s| ^^ + °^^^^A'ii(l-A)<^2 r(di+d2 + l) 
Further, assume for all large a; that 

H{x) = (1 + o(l))Ma;^exp(-r.T*), M > 0,r > 0,6 > 0, N eM. (4.15) 

Since, for any t GlR, 

H{x + tx^~0/{re)) 
^™ WT^ = exp(-t) 



a; — >oo 



we have i7 G MDA{A, w) with 

^(x) = rex^-^, x>0. (4.16) 
In view of Theorem 17.41 the distribution function F of BY satisfies F E AIDA{A,w) and, as a; ^ 00, 

F{x) = (l + o(l))C*a;^-''('^i+''^)exp(-rx''), 

with 

c* = M{re)-^^^-^^^j,^^^nd, + md2 + i). 
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4.2 Regularly Varying Tails 

We deal next with distribution functions p in cither the Frechet or the WeibuU max-domains of attraction. 
As we will discuss below, the asymptotics of i?Q,;3 when H is attracted to the Frechet distribution is quite well 
known, and results for the Weibull max-domain of attraction are less complete. In Section 7 we offer simpler 
proofs of these results, and their converses, i.e., when Ha^p belongs to one of these max-domains of attractions, 
then so does H. 

The unit Frechet distribution function with positive index 7 is $7 (a;) := exp(— x""*"), a; > 0. It is well-known 
that a distribution function H with infinite upper cndpoint rn ^ 00 is in the Frechet max-domain of attraction 
(see e.g., Falk et al. (2004), Resnick (1987)) iff H is regularly varying at infinity with index —7 < 0, i.e., 

hm B^^t-\ yteiO,^). (4.17) 

x^co H[X) 

li Ih = and < 7 < 1, then this condition is the criterion that H is attracted to a positive stable law with 
index 7. Breiman (1965, Proposition 3) shows that if this holds, then the distribution function F oi X = BY, 
where the random multiplier B is independent ofY, is also attracted to the same positive stable law provided 
that < 00. (Thus B is not restricted in sign or magnitude.) Specifically, H and F arc tail equivalent, 

i.e., 

F{x) = il+o{l))E{B''}H{x) (x-^oo). (4.18) 

Jcssen and Mikosch (2006, p. 184) observe that Brciman's proof is valid for any positive 7 if S > and 
^I^T+ej <; 00 for some e > 0. Herman (1992, Theorem 12.3.2) proves this tail equivalence for the case 
B = ^1 - B^^p. 

So in particular, we conclude that if > and H^jj, is defined via (|2.3p with Hajj{Q) = 0, then 

HcAx) = {l + o{l))E{Blu}H{x), x-^oo, (4.19) 

and 

^ "-^^ T{a)r{a + (3 + j)- 

The next theorem asserts that this tail equivalence holds also if 7 = 0, and conversely, if 7 > and Ha.p G 
MDA{^-y), then so is H. Brciman's methodology is completely analytical, and in Section 7 we shall give a 
much simpler proof for the case of a general bounded multiplier < B < 1. We indicate too how this can be 
extended to the general result. 

Theorem 4.4. Let H, H^^p, a, j3 > be two distribution functions defined via (|2.3p with H{0) = 0. Then H 
satisfies (|4.17p with some J > 0, iff Ha^p satisfies (j4.17p with the same index 7. Furthermore, for any 7 > 
we have 

hm = — = 7. (4.20) 

Example 2. Theorem 14.41 shows in particular that Pareto tails are preserved under independent beta random 
scaling. 

The unit Weibull distribution function with index 7 > is 4'^(x) exp(— Ix-p), a; < 0. It is well known that if 
H has a finite upper cndpoint (say r/f = 1), then H E MDA{'^^) iff 

^.^H{l__tx)_ ^ ^ ^^^^^ 

xio H{l-x) 



Theorem 12.3.3 in Bcrman (1992) is closely related to the following result, and in Section 7 we prove a general 
theorem which subsumes both direct assertions. 
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Theorem 4.5. Let H, Ha.i3,a, [3 be as in Theorem \4.4\ If I'm = l,-ff(0) = and (|4.2ip holds for some 7 > 0, 
then Ha43& MDA{^ f3+^) and 

Ha.p{l-x) = {l + o{l))Kx^H{l-x), xjO, (4.22) 

with K := r(a + /3)r(7 + l)/(r(a)r(7 + 13+1). 
Furthermore we have 

xhn r(1 — x) , , 

\n-n=^ f = /3 + 7>0. 4.23 

orio H^,0il-x) 

Conversely, if Ha,p G MDA{^p+j),j > 0, then (|4.2ip is satisfied. 

Remark 4.6. (a) If (|2.3p holds with Ba,i3 ^ gamma ( a, then in Lemma 11 of Hashorva et al. (2007) it is 
shown that H satisfies (|4.17|) with some 7 > 0, iff Ha^jj satisfies (|4.17p with the same index 7 (see also Jessen 
and Mikosch (2006)). 

(b) Under the Gumbel or the Weibull max-domain of attraction assumption on H or Ha, [3 by (j4.5p we have 

lim = 0, 

xTrH H{x) 

whereas when H or Ha,i3 are in the Frechet max-domain of attraction the above limit is a positive constant. 



5 Conditional Limiting Results 

Let the bivariate random vector (Oi, O2) be imiformly distributed on the imit circle, R ^ H he independent of 
(Ol, O2), and let {Si, S2) == R{Oi, O2) be the corresponding bivariate (planar) spherical random vector. Finally, 
define the bivariate elliptical random vector 

{U,V)^{Si,pSi + ^l-p^S2), pe{-l,l). (5.1) 

Distributional properties of spherical and elliptical random vectors are studied by many authors, e.g., Cambanis 
et al. (1981), Fang et al. (1990), Kotz et al. (2000) and their references. Referring to Cambanis et al. (1981) we 
have 

Oj = ^^^1-^(^/2,1/2). (5.2) 

Basic asymptotic properties of spherical and elliptical random vectors can be derived utilising (j5.1[) and (|5.2[) . 
One line of enquiry is to determine the asymptotic behaviour of the conditional distribution ofV — pU given an 
event constraining the values of U . For example, in several statistical applications (see Abdous et al. (2005)) 
the approximation of the conditional random variable 

Z* = {V - px)\U > X, xeM 

is of some interest. Since V — pU = ■\/l — p^5'2, the outcome follows directly from Theorem 12.3.3 in Berman 
(1992), i.e., if i7 G MDA{A,w), then 

c{x)Z* S v^l - p'^Z x^rn, (5.3) 

where c{x) := ^Jw{x) / x, x > 0, and Z is a standard Gaussian random variable. Abdous et al. (2005) is an 
independent account. Theorem 15. II below embellishes this outcome. 
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The point-wise conditioned random variable 

Z.,, = {V ~ px)\U ^x, xgM 

is a particular case of the conditional multivariate models introduced by HefFernan and Tawn (2004) for treating 
certain inference problems. They raise the general problem of conditional limit laws when one component of a 
random vector tends to infinity, and they give results for some particular parametric families. It is known that 
(Hashorva (2006), Corollary 3.1) that Zx has the same Gaussian limit law 

c(x)Z, 4 ^l-p^Z, X t rn. (5.4) 

We will prove that if H is absolutely continuous then (|5.4p holds in the stronger sense that the density functions 
converge. We prove in addition that both limit assertions hold assuming that the (marginal) distribution of \U\ 
is attracted to the Gumbel distribution. Finally, Hashorva and Kotz (2009) gives an account of these results 
based on the strong Kotz approximation. 

Theorem 5.1. Let H,{U,V),p E (-1, 1), c(.t), Z^, Z*, a; > be as above with \U\ - G and G(0) = 0. // 
G G M DA{A,w) or H E MDA{K,w), then (a), (|5.3p is satisfied; and (b), (|5.4p is satisfied if, in addition, H 
is absolutely continuous. 

The proof of this theorem rests on a closure lemma for distributions attracted to the Gumbel distribution. 
Lemma 5.2. Let 0<X^F,p>Obea constant, and denote the distribution function of by Fp. Then 
F e MDA{A,w) iff Fp e MDA{A,Wp) where 

Wpix) = p-'^x^^/P'^-'^w (^x^/P^ . 



6 Estimation of Conditional Survivor and Quantile Function 

For i = 1, 2, . . . , let {Ui, Vi) be independent copies of {U, V) as defined in the previous section, and suppose too 
that R ^ H Cz MDA(A, w) with rn ~ oo. We are interested in the conditional survivor function 

*,(y) := P{V>y\U>x}, x,y eM. 

Estimation of the distribution function 1 — 'l'x(y) when x is large is discussed in detail by Abdous et al. (2007). 
As noted there, if x is large there may be insufficient data available for the effective estimation of "^xiv)- Similar 
difficulties apply for estimation of the inverse function (or conditional quantile function), &{x, ■), s G (0, 1), a; G -2? 
of 1 — ^^xi')- The Gaussian approximation implied by Theorem 15. II entails 



sup 

V<£R 



*:,(yV';^7iiKaO - $(j//v/r^p2) ^0, x^oo, (6.1) 



where $ is the standard Gaussian distribution function. 

On this basis, Abdous et al. (2007) propose two estimators of 'i/x- Theorem 15.11 implies that the Gaussian 
approximation in (|6.ip is valid if we assume instead that U ^ G E AIDA{A,w). For estimation purposes this 
fact is crucial because we can estimate w based only on the random sample C/i, . . . , or Vi, . . . ,Vn- 

A non-parametric estimator /5„ of p is given by (see e.g., Li and Peng (2009))) 

Pn := sin(7rf„/2), n > 1, (6.2) 

where f„ is the empirical estimator of Kendall's tau. Now, if Wn{x) is an estimator of the scaling function w{x) 
(for all large x), then by the above approximation we can estimate ^x{y) by 

^nAv) ¥(/i„(y-/5„x)/(l-p2)i/2^^ (g3) 
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where hn{x) := (it;„(x)/a;)^/^, x > 0. An estimator for the quantilc function O is then given by 



e„(x,s) = pnx + ^l- pI<^ \s)/hnix), a: > 0,s e (0, 1), (6.4) 

with <I>^^ the inverse of <&. Both of these estimators are consequences of the Gaussian approximation. However, 
our concern here is with estimation of w. Specifically, we assume that the scaling function w satisfies (j4.7p with 
positive constants r,9 and Li regularly varying with index 9fi,p, G (— oo,0). It follows that (see Abdous et al. 
(2007)) 

G{x) = exp(-rx^(l + L2(x)) (6.5) 

holds for all large x, where L2 is another regularly varying function with index —dfi. This places G in the class 
of WcibuU-tail distributions, and is the so-called Weibull tail-coefficient (see Gardes and Girard (2006), or 
Diebolt et al. (2007)). Canonical examples of WeibuU-tail distributions arc the Gaussian, gamma, and extended 
Weibull distributions. Next, define for i = 1, . . . , 71, 



R 



and write r['^1 < • ■ ■ < Rir}i,k ~ 1,2 for the associated order statistics. Based on Rl^\i < n or R^i^ ,i < n we 



may construct the Gardes-Girard (2006) estimator of 6, 

n 

where 1 < fc„ < n,r„ > 0,n > 1 are constants satisfying 

lim kn = 00, lim — = 0, lim log(T„/fc„) = 1, lim ■\/fc„6(log(n/fc„)) A eJR, 

n — 'oo n — »oo Jl n — *oo n — 'oo 

and the function b (related to Li) is regularly varying with index rj. The scaling coefficient r can be estimated 
by (see Abdous et al. (2007)) 

f'ri' = -Ly- y^/'),, , ,^l,2,n>l, (6.6) 

leading to the following estimators of w, 

wipix) = X >0,j = 1,2,71 >1. (6.7) 

Our suggestion is to estimate w by wi^' , because it is based on independent and identically distributed Ri,i < n. 
This differs from the estimator Wn recommended by Abdous et al. (2007) which is based on the dependent 

f2l f2) 

random variables i?i„ , . . . , Rmi (recall p„ is estimated from {Ui, T^), i > 1). 

A third estimator of w can be easily constructed by considering the sample Vi, . . . ,Vn since by the assumption 
U = V. 

Note in passing that if = 1, then we have the estimator of r (of interest for G in C{r),r > 0) 



1 ^ \og{n/i) 



7 Further Results and Proofs 



We present first some asymptotic results for the Weyl fractional-order integral operator, followed by the proofs 
of all the results in the previous sections. 
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Theorem 7.1. Let H be a univariate distribution function with H{0) = 0, rn G (0, oo], and H G AIDA(A, w). 
If a is real and P > 0, then 

{Jp,P^^H){x) = (l + o(l))(u;(x))~(^-i)x-"i?(.T), xUh, (7.1) 



{IpP^„H){x) = ii±^(J'^,p_if)(x), xUh. (7.2) 

Proof of Theorem 17. II Let Wx be a random variable whose survivor function is 

nnjr ^ H{x + z/w{x)) 
H{x) 

where 

r(x) = [rn — x)w(x) if rn < oo, & = oo if r/f = oo. 
Then (|4.4p is equivalent to the convergence assertion Wx —* W which has the standard exponential distribution. 



Observe now that for x G (0, r//) we may write 



x~'^H(x) r^""^ ( z \^~^ 

^ ^ ' ' ^ (1 + zMx))-"4P{T^, < z}, 



r(/3) V^(^), 

where v[x) = xw(x) and we have used the substitution y = x + z/w{x) for the second equality. Hence 

{w{x))>^-^x'^ 



iJp,P^.H) = {Wf (1 + Wx/v{x))-"} . 



H{x) ' r(/3) 

It follows from (|4.5[) and the moment convergence theorem (Feller (1971, p. 252)) that the expectation converges 
to E{W'^-'^} ^ r(/3). This proves (TTTjl . 

The same manoeuvres yield 

^ ^ " -il0P^a.H){x) = — - / zP-\l + z/v{xr"P{Wx > Z}dz ^ 1, 



using the dominated convergence theorem, and (|7.2p follows. □ 

Theorem 17.11 subsumes and generalizes results in Berman (1992, §12.2) applying to the case rn = oo. To 
align with Berman's notation, we use /3 — 1 to denote his parameter p, and in what follows we assume that 
EiYf^} < oo. 

(i) Propositions 12.2.3 and 4 in Berman (1992) concern distribution functions F having the form 

POO 

F{x) - (1 + o(l))c / {y- xf-^H{y)dy. 

It is easily seen that 

F[x) = (1 + o{l))cV{(i){Jp+i^j„H){x) = (1 + o[l))cT[(i){w[x)r^H{x), 

and this is valid if /3 > 0, which extends the range of parameter in Berman's Proposition 12.2.4. 

(ii) Proposition 12.2.5 in Berman (1992) concerns survivor functions proportional to the order-g stationary 
excess distribution generated by iJ, i.e., 

/•oo 

F(x) = (l + o(l))c / y^-^H{y)dy, 
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where q is real. The integral can be recast as 

/>OC 

q-^ / {y" - x'i) dH{y) = q-^x''H{x)E{{l + W,/v{x)y - 1}, 
from which it foUows that, as x — > oo, 

F(.)^(l + o(l))c.-^fg). 
(iii) The order-g size-biased distribution generated by H induces survivor functions of the form 
F{x) = (1 + o(l))c / y'^dHiy) = (1 + o{l))c{J,,p^H){x) = (1 + o{l))cxm{x). 
It follows from (|4.1ip that each above F e MDA{A, w). 

Note that if g > 0, then the results under (ii) and (iii) arc related via Theorem 14. II because if Yg and Yg denote 
the order-g size-biased and stationary excess versions of Y, then Yq = YqBq i. See Takes (1996, §4) for this 
connection and further generalization involving beta scaling. 

Theorem 7.2. Let H be a univariate distribution function with rn = oo. Assume that H{0) = and (|4.17p 

holds with 7 > 0. // /3 > and c are two constants such that /? + c < 7 + 1, then 

iJ,,,H){x) ^ (l + oil)):^^^±l^^^Hix)x^'+^-\ x^oo. (7.3) 



Furthermore i/ 7 > 0, then 



{l0PcH)ix) = ^^00. (7.4) 

1 (7 - c) 



Proof of Theorem [Til] Let have the distribution function max(l - H{xt)/ H{x),0). Then (|4.17p is 
equivalent to: If 7 > 0, then Wx — > W which has the Pareto survivor function tr'^ for t > 1; and if 7 = 0, then 
Wx^oo. 

Substituting y = tx into the integral defining J^p.p^H gives the representation 
If 7 > and e > is chosen so/3 + c + e<7+l, then £;(M^'3+c+e-i) 

< 00, and hence the above expectation 

converges to 

E [{W ~ if-^W] = 7^(7 + 1 - /? - c, /3), 

and (|7.3p follows. This assertion follows too if 7 = because {Wx — 1)^~^W^ < VKf+'^~^, and the exponent is 
negative. 

The same substitution yields 

and it is clear that the integral converges to 

{t - if-h'^-'dt = B{j-l3-c, f3), 
whence (Ol). □ 
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Theorem 7.3. Let H be a univariate distribution function with upper endpoint rjy = 1. Assume that H{0) — 0, 
and that (j4.2ip holds with 7 > 0. If (3 > and c€lR are constants and 7 > 0, then 



and if J > 0, then 



{j0,pH){l~x) = (l + o(l))^|^±^i/(l-x)x^-\ xiO (7.5) 



iIpPcH){l~x) = il + o{l))^-x{Jfi,p^H){l-x), xlO. (7.6) 

7 r P 



Proof of Theorem I7.3I Let < 1 be a random variable having the distribution function H{l — tx)/H{l — x). 
If 7 > 0, then (|4.2ip is equivalent to W := U^^'^, where U has the standard uniform distribution (i.e. 

beta(l, 1)), and if 7 = 0, then 1. The substitution y — I xt yields 

{Jp^,H)il -x)^ E]—^^E{{1 - W^r-^l - xW^Y}. 

If 7 > 0, then the expectation converges as x | to 

£;{(l-I^>^-i} =7S(7,/3), 
and if 7 = then it converges to unity. So (|7.5p follows in both cases. 
The same substitution yields 

(/^p,i?)(i -x) = ^^^^^^^""^ J\i - tf-\i xtypiw, < t}dt, 

and the integral converges to ,6(7 + l,/3). Thus 

il,3pjl)il -x) = il + o(l))i?(l - x)x^^,^^^L±}L.^ 
and (TTHl) follows. □ 



Proof of Lemma[57T] Since the first two statements are borrowed from Lemma 2.2 in Pakes and Navarro (2007) 
we show next only statement Hi). Let Y ^ Ba,i3, Ba.\ and -Bq+a,;3-a be independent random variables. For 
any A G (0,/?) we have the stochastic representation (see P-ll) ) 



with Y* ^ Ha.\ another random variable independent of -Bq+a,/3-a. Applying ()2.4p we obtain for any x € (0, rn) 

H^A^) = £(^±^a;"(/^p_„_^i7)(x) 
r(a) 

= H^^^"+'U/3-AP-o-/3i?.,A)(x) 

1 (a + A) 



and 



Consequently 



HcA^) = ^%±^x"(/AP-„-A^f)(.T). 

r(a) 



{Il3P^a-fjH){x) = X^{Ifi^xP-l3{I\P-a-\H)){x), 

and the result follows. □ 
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Proof of Theorem 12.21 The proof follows immediately from Theorem 2.2 in Pakes and Navarro (2007) and 
the identity 



(7.7) 
□ 



Proof of Theorem [nH] The identity (|7.7|) impUes that 

r(a + /3)^„ 



1 = 



r(a) 



x°'{If3P-a-p)ix), Va; e (0,ri/), 



and hence ([^ 



-a; (//3P-a-/3)(2:) ——X°'{If3P-a-pH){x) 



r(a) 
r(a) 



r(a) 

{IpP-a-l3){x) - {Ii3P-a-pH){x) 



= r(« + /?^ ^"(/^p_„_^g)(x), (x-e(0,r^)) 
r(a) 

thus the first result follows utilising further (12. 8|) which holds if H replaces H. 

We show next the second claim. Since H{{]) = 0, Lemma 2.1 in Pakes and Navarro (2007) shows that 

i?a,A(0) = i?a,/3(0) = 0- 



Furthermore, both Ha x and ^ are absolutely continuous and 



X = Y*B. 



a+\,l3-Xi 



with r* - ij„ 



X ^ H, 



a,l3- 



Therefore, in order to show the proof we need to check the assumptions of Theorem 12.21 In our case n = 1, 
hence the condition -ff^"^ = -ff^"^ ~ Ha^p is absolutely continuous is satisfied. Since if^^^ is a density function 
and 5 € [0, 1), then clearly i?^^^ G Further we have H^^^ = i?Q,;3 G X5_q_a-i since i?a,/3 is bounded 

by 1. Applying Theorem 12. 21 for any a; G (0, r//) we may write 

Ha,\{x) = 

L yix -r iJ ) 

r(rv 4- X^l „ r 

(a + \){IsP-a-\-lH a,fj){x) + {Js,p-^-xHa,p){x) 



r!" t m ^°+^(^^^(P-c.-AH../3))(a;) 
1 (a + p) 

r(a + A) ^^^^ 



r(a + /3) 
and the result follows. 

Proof of Corollary 13.21 Letting A ^ in (|3.2p we obtain (recall /q/) := ft.) 



□ 



Consequently, we have 



and in view of 



Ha,p{x) 



-ha,l3{x) 



^^^^z"(/^p_„_^i/)(x), Vx G (0, rn). 
r{a + (3) 



r(a) 

r(a + /3) 

r(«) 



(i?(pa V_„_^i?))(x-), Va- G (0, rn) 



{D{p^If3P-a-pH))ix), yx G (0, rn). 
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Since ha^f} is given by (see (22) in Hashorva et al. (2007)) 

r(a) 

the result follows. □ 



haA^) = -^W7r^^''-\J(3,p.^^,^,H)ix), Vx€(0,r^), (7.9) 



Proof of Theorem 13.31 Let Ba+p^^Xi ~ beta(a + (3i, = 0, . . . , k he independent beta random variables 
independent of X and Y. By the assumptions we may write 

X = YBa,p„ 

= Ba,f3iYBa+l3i,l3o-l3i 

= YiBajj^, with Yi = YaBa+[3^,Po-[3^ = YoBa+pi,Xi, Yq := Y. 



Similarly 



and repeating we arrive at 



X = Y2Ba,p2i with y2 = >lBa+/32_A2 



X = YkBa.f}^, wMhYk ^Yk-iBa+p^^Xk 



Setting Yk+i := X we may write the above stochastic representation as 

Yk+l = YkBa+fJ^^^^Xk + i- 

Let Hq := H and Hk+i := Ha. p. Applying (|3.3p we obtain for any i = 1, . . . , fc + 1. 

H,-l{^) = ^(I^+^i,) ^^"^^-^ [(« + mIs^P~c.^P.~lH^){x) + (j5.,p__,,F,;)(x)] , (7.10) 

and the assertion follows. □ 

We precede our account of scaling relations for the Gumbel distribution with the following proof. 
Proof of Lemma [4?2] If /? > it follows from Theorem mi that 



HiA^) = iJ0+i,p,K^) = (1 + o(i))r(i + T th). 

On the other hand, if B := Bi^^j and c > 1, then 

HiA^) > / P{B > x/y}dH(y) > P{B > c-'^}H{cx). 
Combining these estimates yields 

limsup(xw(a;))'^=- — < oo. 

x^oo H{x) 

The assertion (|4.13p follows by choosing (3 > ji and appealing to (|4.5p in the case rn = oo. □ 



The next result is the foreshadowed generalization of the direct assertion of Theorem 14.11 It comprises two 
parts which respectively yields a tail estimate of the distribution function of a random scaling, and its density 
function. 

Theorem 7.4. Suppose H G AIDA{A,w). (a) If (j){u) > is defined and bounded on [0,1] and it satisfies 

q^{u) = {l + oil))C{l-uf, un, (7.11) 

where /3, C > are constants, then 

I{x) := r c^{x/y)dH{y) - (1 + o(l))CT(l + /?)-ffl-, x T r^. 

Jx [XW[x))t^ 
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(b) If g{u) > is defined in [0, 1] such that ug{u) is defined and bounded on [0,u'] for any u' < 1, and 

giu) = {1 + o{l))c{l - uf-' uU, (7.12) 
where c > and /3 > are constants, then 

J{x) / y-^g{x/y)dH{y) = (1 + o{\))cV{(3) > x ] m. 

Proof of Theorem 17.41 We prove only (b) since the details for (a) are similar and simpler. If v! £ (0, 1), 
then x/y < u' ii y > x/u' and 

/■OO 

Ji{x) := / y-^g{x/y)dH{y) = 0[x~^H{x/u')]. 

J x/u' 

If rn is finite, then Ji{x) = if x > u'rn- If — oo, then, recalling that v{x) = xw{x), Lemma 14.21 ensures 
that Ji{x) ~ o {x^^ Il(x)[v(x))^ (x oo) for all positive /i. 

If c is positive and < e ^ c, then it follows from (|7.12p that g(u)/(l — uY^^ e (c — e, c + e) if u' < u < 1 and 
u' is sufficiently close to unity. Hence J(x) — J\(x) is asymptotically equal to 

Ux) :=c y-\\-xlyf-^dU{y). 



Proceeding as in the proof of Theorem 17.11 we obtain the representation 

The expectation converges to E\W^~^^ = r(/3). Taking [i> (i above, we see that J\(x) = o(J2(a;)), and the 
assertion follows. □ 

Proof of Theorem UTT] Assume that H e MDA{k, w). The direct assertion (|4.9p follow from Theorem Fm a) 
by setting 'f>(u) :=- P{Ba^fi > u} and checking that (|7.1ip holds with C = [PB{a, /3)]^^ . Next, taking g{u) as 
the density function of Ba.fj it is obvious that the conditions of Theorem l7.4r b) are satisfied with c = l/B{a, f3). 
Thus (prU)) follows from gH) and HjTJ^ . 

To prove the converse, assume that Ha^/s E M DA{A, w) for some positive scaling function w. With the notation 
of Theorem 13.31 we may write for i ~ 1 , . . . , /c + 1 



= + A)fep^a-ft-iH»)(x) + (Jg.,p_„^g gO(x)j , yxe{0,rHl7.U) 

where Hq := H, Hk+i H^.p- In view of Theorem 17.11 and (|4.5p . we obtain for i = fc + 1 that 

H,.^{x) = (l + o(l))-I^^-t^x"+ft-(J,,,,_„_,_iJ,)(x) 

= {l + o{l))^^^^p^,x^-^-^^{w{x)r^'^-m{xl Vx-Tr«. 

By ((43)) and (|46)) it follows that Hk G A/i:)A(A, w). Since the above holds for alH = 1, . . . , fc, it follows that 
Ho = H e MDA{A,w) too. Next, jLH) and jLH) imply for any x > that 

ha,(3{x) ^ x{Jf)^p_^_^,H){x) 

Ha,0{x) {Ipp^c-p-iH){xy 
so applying Theorem 17. II establishes (j4.10p . and the result follows. □ 
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As foreshadowed above, the following argument includes a simple proof of (|4.18p for an arbitrary bounded 
random scaling. We then show how this proof can be extended to remove the boundedness restriction. 
Proof of Theorem 14.41 With Wx as in the proof of Theorem 17. 2[ clearly 

P{XB >x} = H{x)E{P{B > W^^}}. 

But ^ 

P{B > W.^^} P{B > W-^} / P{W > u^^}dP{B <u} = E{B^}, 

Jo 

and hence F{x) = (1 + o{l))H{x)E{B'^}. 

Similarly, if B has the density function g{u), then the density function of F is 

fix) = x-m{x)E{W-'g{W-')}. 

If g satisfies appropriate boundedness conditions, which certainly are satisfied by beta density functions, then 
the expectation converges to 



E{W''g{W-')} 



hjiw-'^)dw =-fE{B'^}. 



It follows that {h{x)/H{x)) = (1 + o{\)){'^ / x). The direct assertions of Theorem 14.41 follow. 

The converse asserts that if Hap is regularly varying with index —7 < 0, then H is also regularly varying with 
index —7. If 7 > 0, then the proof follows from Theorem 13.31 and Theorem 17.21 

Alternatively, write can write Ba^p = Z\I(Z\ + Z'l). where y, Zi,Z2, are independent random variables such 
that Z\ ~ gamma(a, 1) and Z2 ^ gamma(/3, 1). Since Z\ ^ Z2 ^ gamma(Q: + /3) is independent of B^jj, the 
relation X = YBa.p is equivalent to X{Zi + Z2) = YZi. It follows from Jessen and Mikosch (2006, Lemma 
4.2(a)) that the survivor function of Y Zi is regularly varying with index —7, and Lemma 17 in Hashorva et al. 
(2007) implies the same is true for H{x). We emphasize that this proof is valid for 7 > 0. □ 

Note that Theorem 12.3.2 in Berman (1992) follows from the above direct proof since 

E{il Ba.p) } - E{B^, J _ 

The situation where B is allowed to be unbounded can be handled by writing 

P{YB > x} = P{YB > x-Y > x} + P{YB > x;Y < x}. (7.15) 

Exactly as in the last proof, the first term on the right is asymptotically proportional to H{x)P{B > ly^^}, 
but now the probability term evaluates as 

P{W > B^^} = P{B > 1} + E{B'^] B < 1}. 

The second term on the right-hand side of (|7.15p equals 

/oo 
(H{x/z) - H{x)) dP{B < z} 

= (1 + o{l))H{x) [E{B^; B>l}-P{B> 1}] , 

provided the limit here can be taken inside the integral. This is permissible if E{B'^^^} < 00 for some e > 0. 

The converse tail equivalence statement is open in general, but see Hashorva et al. (2007) for the case where B 
has a gamma distribution. 

The following result is the analogue of Theorem 17.41 for H E MDA{'i>~f) and it generalizes the direct assertion 
of Theorem HSl 
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Theorem 7.5. Let H{0) = 0, rn ^ 1 and H e ]\IDA{^^). (a) If (f7UT|) holds, then 

iix) = J^' ^ ^{x/y)dH{y) = (1 + o{l))C^^^^^^^^^^x^H{l - x), (x i 0). 



(b) If ({71^ /lo/ds, th 



en 



J{x)=f y-i5(x/y)diJ(y) = (l + o(l))c^Ml±i)^/5-i7J(i_,,), (^10). 
Ji-x -l-(P + 7j 

Proof of Theorem 17.51 For (a) simply observe that if 1 — x < y < 1, then 

1-^) =(l + o(l))C(y-x)V^ (xiO). 
V J 

Hence /(x) is asymptoticahy equal to CT(/3 + l)(i/;3+i,p_^-ff)(l — x), and the assertion follows from Theorem 
17.31 Similarly, j(x) is asymptotically equal to cr(/3)(j7/3,p_^iJ)(l — x). □ 

Proof of TheoremIM] If i7 e MDA{;^^), then ([42211 and (j4:23l) follow from TheoremO 

Proof of Lemma [O] It follows from (|4.5p that 

(y + = (1 + o{l))pyP-'/w{y) = (1 + o(l)) (y^ + {t/wpiyP))) , 



and hence that the necessary and sufficient condition (|4.4[) applied to F is equivalent to 

lim P{XP >yP + t/wpiyP)\XP > yP} = e"*. 

>CX3 

Setting X ^ yP shows this is equivalent to Fp E MDA{A, Wp). □ 
Proof of Theorem 15.11 Let G2 and II2 denote the distribution functions of and i?^, respectively. It 



follows from (|5.2|) . Lemma [5.21 and Theorem 14.11 that 

H e MDA{A, w) iff H2 e MDA{A, W2) iff G2 e A/i:>A(A, wa) iff G G ML>A(A, lo). 

This, together with Theorem 12.3.3 in Berman (1992) implies that ([O]) holds if G G MDA{A,'w), i.e. (a) is 
valid. 

By the same reasoning, (b) follows if we prove it assuming H E MDA{A, w) and H has a density function h. 
Observing that has the same distribution as ■\/ 1 — 5*2 1 Si ~ x, we set p = without loss of generality. In 
this case (following the example of Abdous et al. (2005)) we can use the equivalent representation ([/, V) = 
(/i Vb, /2 V ( 1 — B'^), where Ii, I2 and B are independent, the Ij = ±1 with equal probability, and B ^ 
beta(l/2, 1/2). The joint density function f{u,v) of (U,v) is radially symmetric and a routine computation 
yields 

f{u,v) ~ 2tt\/ V? + h + . 

It is more expedient to work directly in terms of ^12(2^) = (2%/^) ^ ^ (v^); whence 

Integration with respect to v and using the substitution y = v'^ gives the marginal density function of U, 

1 r°° 

fu{u) = - h2{y + u^)y-'/^dy, 



and hence the density function of is 

h2 (x2 + ii2) 



fiv\x) := f{x,v)/fu{x) 



S^h2{y + v'^)y '^I'^dy' 
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valid for real v and x > 0. Note that the distribution of is symmetric about zero. 

Let t > and replace v with t/c{x) in this density function. Since c^(x) = 2w2(x'^), the density fimction of 
c{x)Zx is the function of s = given by 

Divide the numerator and denominator of the right-hand side by W2{s)H2{s). Since /i2(s) = {l+o{l))w2{s) H 2(3) , 
it follows from Lemma [5.21 and (|4.4p and (|4.6p applied to H2, that the numerator term obtained from (|7.16p 
converges to exp(— as x — > r/j. 

Next, making the substitution z ~ yw2{s) in the integral at (|7.16p . the denominator term obtained from the 
division operation is 

^ 

W2{s)H2is) 

where Ws is as defined in the proof of Theorem 17.11 (s replacing x there) . The moment convergence theorem 
ensures that 



^lim E {W-^/^} = E = r(l/2) = V^. 



Combining these limits shows that CiA^) converges to the standard Gaussian density function, and the assertion 
follows. □ 
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